Mon. Not. R. Astron. Soc. 000. HIO (20051 Printed 4 February 2008 (MN WF&t style file v2.2) 



The formation history of elliptical galaxies 

Gabriella De Lucia 1 *, Volker Springel 1 , Simon D. M. White 1 , Darren Croton 1 ! 
and Guinevere Kauffmann 1 

1 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, D-85748 Garching, Germany 



4 February 2008 



in 
o 
o 

(N 

> 
O 

(N 

(N 
> 

in 

(N 

On 
O 

in 
o 

Oh 

6 
& 



ABSTRACT 

We take advantage of the largest high-resolution simulation of cosmic structure growth 
ever carried out - the Millennium Simulation of the concordance ACDM cosmogony 
- to study how the star formation histories, ages and metallicities of elliptical galaxies 
depend on environment and on stellar mass. We concentrate on a galaxy formation 
model which is tuned to fit the joint luminosity/colour/morphology distribution of 
low redshift galaxies. Massive ellipticals in this model have higher metal abundances, 
older luminosity-weighted ages, shorter star formation timescales, but lower assem- 
bly redshifts than less massive systems. Within clusters the typical masses, ages and 
metal abundances of ellipticals are predicted to decrease, on average, with increasing 
distance from the cluster centre. We also quantify the effective number of progenitors 
of ellipticals as a function of present stellar mass, finding typical numbers below 2 for 
ill* < 10 11 Mq, rising to ~ 5 for the most massive systems. These findings are consis- 
tent with recent observational results that suggest "down-sizing" or "anti-hierarchical" 
behaviour for the star formation history of the elliptical galaxy population, despite the 
fact that our model includes all the standard elements of hierarchical galaxy formation 
and is implemented on the standard, ACDM cosmogony. 

Key words: galaxies: formation - galaxies: evolution - galaxies: elliptical and lentic- 
ular, cD - galaxies: bulges - galaxies: stellar content 



1 INTRODUCTION 

Elliptical galaxies are the most massive stellar systems in 
the local Universe and appear to define a homogeneous class 
of objects with uniformly old and red populations, negligi- 
ble amounts of gas, and very little star formation. Their 
deceptively simple appearance inspired a 'classical' forma- 
tion scenario in which they form in a single intense burst 
of star formation at high redshifts (z > 5), followed by 
pass ive evolution of their stellar populations to the present 
day dPartridge fc Peeblesll967l:lLarsorJll975h . This so-called 
monolithic scenario successfully explains the tightness of the 
fundamental scaling relations that elliptical galaxies obey, 
like the colour-magnitude relation and the fundamental 
plane, as wel l as the evolution of t h ese relations as a function 
of re dshift (iKodama et alJ Il998t Ivan Dokkum fc Stanford! 
120031) . 

This classical view has been tenaciously resistant to 
challenges by other theoretical models, despite numer- 
ous indications for a more complex formation scenario. 
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iToomre fc Toomrd il972f) suggested that elliptical galax- 
ies can form from major mergers of massive disk galaxies . 
Detailed numerical simulations jFarouki fc Shapiro! fl983 : 
iNeeroponte fc White! Il983h later showed that the merger 
of two spiral galaxies of comparable mass can indeed pro- 
duce a remnant with structural and photometric proper- 
ties resembling those of elliptical galaxies. In more recent 
years, a large body of observational evidence has been 
collected that demonstrates that interactions and merg- 
ers indeed represent a common phenomenon at high red- 
shifts, and that these processes affect the population of el- 
liptica l galaxies in the local Universe. ISchweizer fc Seitzer! 
( 1992) found evidence for bluer colours of elliptical galax- 
ies with increasing morphological disturbance in a study 
based on a small sampl e with a strong bias towards iso- 
lated systems (see iMichard fc Pru gnicl 2004$) • Later stud- 
ies using absorption-line indices have demonstrated that a 
significant fraction of cluster early-type gal axies has un- 
dergo ne recent episodes of star formation jBarger et alJ 
1996). Signs of recent star formation activity have also 
been detected in a n umber of high redshift early-type 
galaxies using colours jMenanteau, A braham fc Ellis! l200ll : 
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and emission line diagnostics jTreu et alJl2002l : IWillis et alJ 
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2002). These results favour, at least for a part of the ellip- 
tical galaxy population, a hierarchical formation scenario in 
which larger spheroidals are assembled relatively late from 
the merger of late-type galaxies of comparable mass. Such 
a bottom-up formation scenario is naturally expected for 
the structure formation process in cosmologies dominated 
by cold dark matter. 

However, despite an enormous amount of work both on 
the theoretical and on the observational side, the debate 
concerning the two competing theories for the formation of 
elliptical galaxies has remained open. As described above, 
a relatively large fraction of early-type systems shows clear 
evidence of interactions, mergers, and recent star formation. 
However, the data also seem to indicate that only a small 
fraction of the mass is involved in such episodes. The lat- 
ter observational result has often been interpreted as strong 
evidence against the more extended star formation history 
naively predicted from hierarchical models. A related issue 
concerns the a-element enhancements observed in ellipticals. 
The so-called a-elements are released mainly by supernovae 
type-II, while the main contribution to the Fe-peak elements 
comes from supernovae type-la. For these reasons, the [a/Fe] 
ratio is believed to encode important information on the 
time-scale of star formation. It is now a well established re- 
sult that massive ellipticals have super-solar [a/Fe] ratios, 
suggesting that they formed on relatively short time-scales 
and/or have an initial mass function that is skewed towards 
massive stars. The inability of early models of the hierar- 
chical merger paradigm to reproduce this observed trend 
has been point ed out as a serious problem for these models 
llThomaslll999h . 

Another contentious issue is related to the significantly 
dissimilar evolution of early-type galaxies in different envi- 
ronments p redicted by early semi-analytic models of galax y 
formation jKauffmannlll996bl : iBaueh. Cole fc Frenklll99fi) . 
Such differential evolution is a natural outcome of the hi- 
erarchical scenario, because present day clusters of galaxies 
form from the highest peaks in the primordial density fields, 
leading to an earlier onset of the co llapse of the dark matter 
haloes and to more rapid mergers jKauffmannll 19951) . 

In recent years, considerable observational progress has 
been made in the study of the stellar populations of el- 
liptical galaxies in dif ferent environments and at differ- 
ent redshifts fe.e. iThomas et aljEfjoS iDenicolo et alJl2005l; 
iTreu et aljf2005l: Ivan der Wei et al.ll2005h . One has to bear 
in mind however that the derived ages and metallicities de- 
pend quite strongly on the models employed in the analy- 
sis, primarily because of an essentially unavoidable intrinsic 
age-metallicity degeneracy, which ap pears to be stronger fo r 
older and more metal rich systems iDenicolo et al.l 12005). 
The situation is especially unclear for 'field' galaxies, since 
here the very definition of 'field' is fraught with ambiguities. 
Significant systematic uncertainties are also present in theo- 
retical models used for studying galaxy formation, where 
poorly understood physical processes need to be treated 
with coarse approximations in order to predict observable 
properties of galaxies. 

A detailed analysis of the properties of elliptical galaxies 
expected in the framework of hierarchic al galaxy formation 
has been given in a number of papers (|KauffmanrJ ll996b: 
iBaugh et af]ll99fi iKauffmann fc Charlodll99^ . This early 
work used an extension of the Press-Schechter theory to pro- 



duce Monte Carlo realizations of the merger trees of dark 
matter haloes, thus allowing the progenitors of a dark matter 
halo to be followed back in time to arbitrarily high redshifts. 
A drawback of this approach is its lack of spatial informa- 
tion on the clustering of galaxies. Another lies in the inher- 
ent inaccuracies of the Press-Schechter formalism. Recent 
years have witnessed substantial progress in this regard with 
the advent of "hybrid" techniques where high-resolution 
N-body simulations of structure formation are used to di- 
rectly measure dark matter merger history trees from sim- 
ulations which are then combined with s emi-analytic sim- 
ulations of the galaxy formation physics (IKauffmann et alJ 
119991: ISprineel et aLlboOlt iMathis et al.ll2002l) . This allows 
the spatial and kinematic distribution of model galaxies to 
be predicted as a function of redshift, providing a more di- 
rect and more powerful comparison between theoretical pre- 
dictions and observational results. Also, this approach elim- 
inates much of the uncertainty introduced by Monte-Carlo 
prescriptions for producing mock merging histories. 

Most previous semi-analytic studies of the properties 
of ellipticals were carried out in the framework of a cosmo- 
logical model with critical matter density. This cosmogony 
has been replaced in recent years by the ACDM scenario, 
which has become the de facto standard cosmological model, 
thanks to its concordance with a variety of observational 
data, including the most recent cosmic microwave back- 
ground measurements, distant supernova observations, and 
cosmic shear measurements. Given that both the cosmolog- 
ical model and the modelling techniques have changed sig- 
nificantly, it is interesting to revisit the question of the for- 
mation and evolution of elliptical galaxies. In this paper, we 
present the results of applying a semi-analytic model that 
tracks dark matter halos and their embedded substructures 
to the largest high-resolution simulation of cosmic structure 
growth ever carried out. Here we concentrate on the analysis 
of the star formation histories, the ages, and the metallici- 
ties of model elliptical galaxies as a function of galaxy mass 
and of environment. In a companion paper, we will study 
how the distribution of metallicities and ages depends on 
the feedback model and on the chemical enrichment scheme 
employed, exploring also the relative contribution of these 
two in shaping observed scaling properties like the colour- 
magnitude relation. 

The layout of this paper is as follows. In Section [5] we 
briefly describe the simulation used in our study, and in Sec- 
tion|^]we give a concise overview of the semi-analytic model 
employed for our analysis. In Section [I] we discuss how the 
star formation history of model elliptical galaxies depends 
on the stellar mass of galaxies and on the environment, while 
Section [S] discusses the dependence of ages and metallicities 
on galaxy mass and environment. Finally, in Section we 
summarise and discuss our findings, and give our conclu- 
sions. 



2 THE SIMULATION 

In this study, we analyse a large collisionless cosmologi- 
cal simulation which follows N — 2160 3 particles of mass 
8.6 x 10 8 ^'Mq within a comov ing box of size 500 /i _1 Mpc 
on a side ijSpringel et al.l l2005Tl . The spatial resolution is 
5ft _1 kpc, available everywhere in the periodic box. The 
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cosmological model is a ACDM model with parameters 
Sl m = 0.25, Q b = 0.045, h = 0.73, Q A = 0.75, n = 1, and 
as = 0.9, where the Hubble constant is parameterised as 
Ho = 100 /ikms _1 Mpc _1 . These cosmological parameters 
are consistent with recent determinations from the combined 
analysis of the 2dFGRS (Collcss ct al. 2001) and first year 
WMAP data JSpergel et al.ll2003ft . 

A remarkable aspect of this Millennium Simulation car- 
ried out by the Virgo Consortium 1 is its good mass reso- 
lution combined with a very large particle number, more 
than 10 billion. It is the largest high-resolution simulation 
of cosmic structure growth carried out so far. This provides 
substantial statistical power, sampling the format ion history 
even o f rare objects in a representative fashion. iGao et all 
(2005) have exploited this to show that the clustering of ha- 
los of fixed mass does depend on their formation redshift, 
an effect that is weak for massive systems but becomes pro- 
gressively stronger towards smaller galaxies. We note that 
simpler schemes for modelling the galaxy distribution, based 
for example on the halo occupation distribution schemes or 
on Monte Carlo halo merger trees, do not account for this 
effect, highlighting the importance of direct simulations of 
structure formation for obtaining fully reliable merger his- 
tories. 

The simulation has sufficient resolution to track the mo- 
tion of dark matter substructures in massive halos, mak- 
ing it possible to follow the orbits of cluster ellipticals in 
an unambiguous fashion. During the simulation, 64 time 
slices were saved, together with group catalogues and their 
embedded substructures, the latter defined as locally over- 
dense and self-bo und structures, iden tified with the SUB- 
FIND algorithm dSpringel et al.ll200ll) . These group cata- 
logues were then used to construct detailed merging history 
trees of all gravitatio nally self-bound dark matter structures 
fepringel et al.ll2005l) . The merger trees describe the assem- 
bly of about 20 million galaxies, and form the basic input 
needed by the semi-analytic simulation of the galaxy forma- 
tion process considered in this study. 



3 THE SEMI-ANALYTIC MODEL 

Our technique for grafting the semi-analytic model onto the 
Mille nnium Simulati on is similar in spirit to that used by 
ISprineel et all (1200 ll) and De Lucia, Kauffmann & White 
(2004), but has been updated in a number of impor- 
tant points. A full descr iption of the mode l is given by 
ISprineel etHI <l2005l) and lCroton et alJ J2005I) . but we here 
give a brief account of those aspects that are particularly 
relevant for the present study. 

One of the key differences of our approach from tra- 
ditional semi-analytic models is that we explicitly follow 
dark matter haloes even after they are accreted onto larger 
systems. This allows the dynamics of satellite galaxies resid- 
ing in the infalling haloes to be properly followed until the 
parent dark matter 'substructure' is compl etely destroyed 
because of tidal truncation and stripping llDe Lucia et alJ 



1 The Virgo consortium Ihttp:/ /www. virgo.dur.ac.uk/ I is an in- 
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l2004al : lGao et al.ll2004tl . At this point, we estimate a resid- 
ual survival time of the satellite galaxy and track its po- 
sition by means of the most bound particle of the subhalo 
identified just before the substructure was disrupted. In this 
work, we consider as genuine substructures those subhalos 
that contain at least 20 self-bound particles. The parent 
catalogue of dark matter haloes, which are analysed for sub- 
structures and decomposed accordingly, is identified with a 
standard friends-of-friends (FOF) algorithm with a link- 
ing length of 0.2 in units of the mean particle separation. 
In our model, we assume that only the galaxy located at 
the position of the most bound particle of the FOF halo - 
the 'central galaxy' - is fed by radiative cooling from the 
surrounding halo, i.e. genuine satellite galaxies cannot re- 
plenish their reservoir of cold gas. Gas infall and coolin g are 
modelled as described in detail in lCroton et alJ J2005). 

Hierarchical merging trees form the backbone of our 
semi-analytic model, and are built for all the self-bound 
dark matter halos and subhalos in the Millennium Sirmi - 
lation using the methods described in lSpringel et alJ i2005l) : 
the descendant in the next time slice of each dark matter 
(sub)halo is identified as the (sub)halo that contains the 
largest number of its most tightly bound particles. Individ- 
ual trees are stored separately in a self-contained fashion, so 
that the semi-analytic code can be run for each of these trees 
sequentially, instead of having to process the whole galaxy 
population in one single run for the entire simulation box. 
This makes the computation feasible even on small work- 
stations and allows for easy parallelisation, such that the 
computation of the galaxy properties can be repeated with 
different physical assumptions in a matter of hours, if de- 
sired. 

The descriptions we adopted for modelling the various 
mixing and exchange processes occurring between different 
galactic phases (stars and cold gas, hot diffuse gas in the 
dark matter haloes, and intergalactic gas outside virializcd 
haloes) are essentially the same as those in iDe Lucia et all 
(2004bJ). In the present study, we use their 'ejection slow' 
feedback model, which has been shown to reproduce both 
the observed relation between stellar mass and cold phase 
metallicity, and the relation between luminosity and cold gas 
fraction for galaxies in the local Universe, as well as the ob- 
served decline in baryon fraction from rich clusters to galaxy 
groups. As in IDe Lucia et~ai] d2004rl . we use metallicity- 
dependent cooling rates and luminosities, and we refer the 
reader to the original paper for more details on these imple- 
mentations. 

In the present paper, we adopt new parameterisations 
of star formation and of the suppression of co oling flows by 
centra l galaxy AGN activity, as introduced bv lCroton et alJ 
( 2005). In the following, we briefly summarise the main char- 
acteristics of these physical prescriptions where they differ 
from those us ed in our previous w ork. 

Following iKauffmannl dl996al) and lCroton et alJ <2005ll . 
we assume that the star formation occurs with a rate given 
by: 

1p = a(M coM - M crit )Adyn, (1) 

where M co id and tdyn = i?disc/Kdr are the cold gas mass and 
the dynamical time of the galaxy, respectively. The dimen- 
sionless parameter a regulates the efficiency of the conver- 
sion of gas into stars. Star formation is allowed to occur only 
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if the gas surface density is larger than a critical value (that 
is used to obtain M cr it in Eq. given by: 

E - - " >< 10? (200^ ) (life) _1 MQ kPC " 2 (2) 

Note that in iDe Lucia et"afl j2004bl) we did not assume a 
surface density threshold for the star formation but we as- 
sumed a dependence of the star formation efficiency on the 
circular velocity of the parent galaxy. This assumption was 
responsible for delaying the star formation in small haloes 
until lower redshifts, which correctly reproduces the ob- 
served trend for incr easing gas fraction at lower luminosi- 
ties. As explained in iKauffmannl il996aT) . the introduction 
of a surface density threshold also naturally reproduces the 
observed trend of the gas fraction as a function of galaxy lu- 
minosity, due to the fact that the gas density always remains 
close to the critical gas surface density value. We postpone 
a more detailed investigation on how the census of ages and 
metallicities is influenced by the assumptions on the star for- 
mation and feedback model to a forthcoming paper. With 
respect to the trends presented in this paper, we have veri- 
fied that the two different assumptions produce very similar 
results. 

As in lDe Lucia et alJ l|2004b|) . we assume that bulge for- 
mation takes place during mergers: in the case of a 'minor' 
merger, we transfer the stellar mass of the merged galaxy to 
the bulge of the central galaxy and update the photometric 
properties of this galaxy. The cold gas of the satellite galaxy 
is added to the disk of the central galaxy and a fraction 
of the combined cold gas from both galaxies is turned into 
stars as a result of the merger. Any stars that formed during 
the burst are also added to the disk of the central galaxy. If 
the mass ratio of the merging galaxies is larger than 0.3, we 
assume that we witness a 'major' merger that gives rise to a 
more significant starburst and destroys the disk of the cen- 
tral galaxy completely, producing a purely spheroidal stel- 
lar distribution. Note that the galaxy can grow a new disc 
later on, provided it is fed by an appreciable cooling flow. 
Our starburst implementation is based on the 'collis ional 
starburst' model int roduced bv lSomerville et all i200ll) (see 
ICroton et all ( 20051) for details). 

Following Croton et alJ i2005l) . we extend the spheroid 
formation by assuming that bulges can also grow from disk 
instabil ities, based on the the anal ytic model for disk forma- 
tion by iMo. Mao fc White! Jl998ft . We note that the addi- 
tion of this chanel for bulge formation does not substantially 
modify the trends presented in this paper. However, this ad- 
ditional physical mechanism changes the relative fractions of 
different morphological types. For the model explored in this 
paper, the final fractions of ellipticals, spirals, and lenticu- 
lars brighter than —18 in the V-band are 17, 65, and 18 per 
cent, respectively, and are very close to the obser ved rela- 
tive fr actions 13, 67, and 20 per cent measured bv lLovedavl 
( 1996). If bulge growth through disk instabilities is switched 
off, the above fractions become: 7, 84, and 8 per cent re- 
spectively. Considering only galaxies brighter than —20, the 
fractions cited above become: 23, 58, and 19 for our default 
model and 13, 67, and 20 for a model where bulge growth 
through accretion is switched off. These numbers suggest 
that this channel of bulge formation may be more impor- 
tant for fainter ellipticals. We will come back to this issue 
in Sec.|S] 



As in previous work, we determine the morphol- 
ogy of our model galaxies by using the B-band bulge- 
to-disc ratio together w i th the observational relation by 
ISimien fc de Vaucouleursl lll986t) between this quantity and 
the galaxy morphological type. For the numbers quoted 
above and for the following analysis, we classify as ellip- 
ticals all galaxies with AM < 0.4 (AM = Mb u i ge — Mtotal), 
as spirals or irregulars all galaxies with AM > 1.56, and 
as lenticulars (SO) all galaxies with intermediate value of 
AM. We include in our analysis all galaxies with stellar 
mass larger than 3 x 10 s Mq. Note that, although we are 
essentially complete down to this mass limit, the less mas- 
sive galaxies included in this analysis reside in substructures 
whose mass accretion history can be followed back in time 
in many cases only for a small number of snapshots. The de- 
termination of the morphological type then becomes quite 
noisy and uncertain. We estimate that our morphological 
type determination is robust for galaxies with mass equal 
to a few times 10 9 Mq . The inclusion of galaxies below this 
limit, however, does not affect our main results of this study. 
When necessary, we will explicitly show results including 
only galaxies with mass larger than 4 x 10 9 Mq , whose mor- 
phological type can be considered 'secure'. Our final sample 
contains 1, 031, 049 elliptical galaxies with stellar mass larger 
than 4 x 10 9 Mq. 810, 486 of these have stellar mass larger 
than 1 x 10 10 M . 

Finally, we use the model of ICroton et alJ |20 05) to de- 
scribe central heating by AGN in massive groups and clus- 
ters and the associated suppression of cooling flows. In this 
model, gas condensation in massive systems is efficiently 
suppressed by 'radio mode' outflows that occur when a mas- 
sive black hole finds itself at the centre of a static hot gas 
halo. The importance of these outflows grows with decreas- 
ing redshift and with the mass of the system. We refer to the 
original paper for full details and the physical motivation for 
this feedback implementation. We will later discuss in more 
detail the effect that this particular implementation has on 
the trends presented in our study. 



4 THE STAR FORMATION HISTORY OF 
ELLIPTICAL GALAXIES 

As discussed earlier, the uniformly red and old stellar pop- 
ulations of elliptical galaxies have traditionally been inter- 
preted as evidence for a formation scenario in which these 
galaxies form in a single intense burst of star formation 
at high redshift and then passively evolve to the present 
day. Direct observations of the formation and evolution of 
early type galaxies are, howeve r, difficult, and are plagued b y 
the so called 'progenitor-bias' llvan Dokkum fc Franxll996ll . 
One approach that is adopted to constrain the formation 
mechanism of these galaxies is that of studying in detail 
their stellar populations by means of population synthesis 
te chniques . This method has its roots in a pioneering study 
bv lTinslevI lll972l) . and has become increasingly popular af- 
ter the introduction of more detailed population synthesis 
mode l s feruzual Aj|l983| : iGuiderdoni fc Rocca-Volmerangd 
Il987t iBuzzonil ^j89fT Recent improvements have come 
through the development of medium to high resolution 
spectral models that include quite complete libraries of 
stellar spectra and improved treatments of stellar evolu- 
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Figure 1. Average star formation histories of model elliptical 
galaxies split into bins of different stellar mass, normalised to the 
total mass of stars formed. The two panels are for galaxies resid- 
ing in haloes of different mass, as indicated by the labels. In both 
panels, the solid line shows the average star formation history for 
all the elliptical galaxies in the sample under investigation. The 
long dashed, dash-dotted, dashed, and dotted lines refer to galax- 
ies with stellar mass ~ 10 12 , 10 11 , 10 10 , and 10 9 M respectively. 
The vertical line in both panels is included to guide the eye. 



tionary theory JVazdekisl l200lt iBruzual fc CharlotJ l2003t 
iThomas. Maraston fc Benderll200:fl. 

The use of these more sophisticated models, together 
with the acquisition of better and larger amounts of 
data have recently established firm evidence for a mass- 
dependent e volutionary his t ory of the elliptical galaxy 
population toe Lucia et all l2004d: iKodama et alJ 12004 
[ Thomas et all 120051 Ivan der Wei et alJ 120051 : iTreu et all 
2005). The data suggest that less massive ellipticals have 
more extended star formation histories than their more mas- 
sive counterparts, giving them a lower characteristic forma- 
tion redshift, in marked contrast to naive expectations based 
on the growth of dark matter halos in hierarchical CDM cos- 
mologies. This observational finding is compatible with the 
"down-sizing" scenario for star form ation proposed earlier 
jFaber et al.lll995l: ICowie et alJll996ft . 

In this section, we study in detail the star formation 
histories of model elliptical galaxies and their dependence 
on stellar mass and environment. In Fig.0 we show the av- 
erage star formation history, normalised to the total mass of 
stars formed, for model elliptical galaxies split into different 
bins of stellar mass. The two panels are for galaxies resid- 
ing in haloes of different mass. In the following we will refer 
to galaxies in haloes with mass > 8 x 10 14 M as 'cluster' 
ellipticals and as 'field' ellipticals to all model ellipticals re- 
siding in less massive haloes. In both panels, the solid line 



shows the average star formation history for all the elliptical 
galaxies in the sample under investigation. The long dashed, 
dash-dotted, dashed, and dotted lines refer to galaxies with 
stellar mass ~ 1 X 10 12 , 1 x 10 11 , 1 x 10 10 , and 1 x 10 9 M , 
respectively 2 . The vertical lines are included to guide the 
eye and mark the peak of the 10 12 MQ-ellipticals in the top 
panel. 



Fig. shows the most important result of this paper: 
more massive elliptical galaxies have star formation histories 
that peak at higher redshifts (~ 5) than lower mass systems, 
and can reach star formation rates up to several thousands 
of solar masses per year for galaxies ending up in overdense 
regions. Less massive elliptical galaxies have star formation 
histories that peak at progressively lower redshifts and are 
extended over a longer time interval. 



A comparison of the top and bottom panels of Fig. 
shows that the qualitative behaviour for 'field' and 'cluster' 
ellipticals is the same, but that for fixed mass, the star for- 
mation histories of field ellipticals are predicted to be more 
extended than those of ellipticals in clusters. This is a nat- 
ural outcome of the hierarchical scenario, where haloes in 
regions of the Universe that are destined to form a cluster 
collapse earlier and merge more rapidly. The star forma- 
tion histories shown in Fig. Q represent averages computed 
over all the elliptical galaxies in the simulation box, but the 
trends remain true also when a much smaller volume of the 
simulation, and hence a much smaller sample size, is anal- 
ysed. Fig. |5| shows the star formation histories of randomly 
selected elliptical galaxies in different mass bins and in dif- 
ferent environments. The figure shows that individual star 
formation histories display a much more 'bursty' behaviour 
than those shown in Fig.Q This reflects our assumption that 
bulge formation takes place during merger-induced bursts, 
which naturally gives the star formation histories of individ- 
ual systems a bursty nature quite different from the smooth 
history seen for the population average. We will comment 
more on the implications of this for the scatter of the ages 
for the model elliptical galaxies in the following section. 



In Fig. |21 we show the star formation histories again, 
but this time split into bins of different parent halo mass. 
The long dashed, dash-dotted, dashed, and dotted lines are 
for elliptical galaxies in haloes of mass ~ 1 x 10 15 , 1 x 10 14 , 
1 x 10 13 , and 1 x 10 12 Mq respectively 3 . Only galaxies with 
stellar mass larger than 4 x 10 9 Mq are used here. The solid 
line shows the average mass-weighted star formation history 
for all the galaxies in the sample. The faster evolution of 
proto-cluster regions produces star formation histories that 
peak at higher redshifts for galaxies in more massive haloes. 
Given that galaxies of a fixed stellar mass occur in haloes 
covering a wide range of masses, it is not surprising that 
the dependence of the star formation history on halo mass 
is much weaker than that on galaxy stellar mass. 
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Figure 2. Each panel shows the star formation history for three randomly selected model elliptical galaxies in the given mass bin. 
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Figure 3. As in Fig.Q but with model elliptical galaxies split into 
bins of different parent halo mass. The solid line shows the average 
mass— weighted star formation history for all the elliptical galaxies 
in the sample under investigation. The long dashed, dash-dotted, 
dashed, and dotted lines refer to galaxies residing in haloes with 
mass M 20 o ^ 10 15 , 10 14 , 10 13 , and 10 12 M , respectively. 



5 THE DISTRIBUTION OF AGES AND 
METALLICITY 

We now turn to an analysis of the distribution of ages and 
metallicities of model elliptical galaxies as a function of stel- 
lar mass and environment. In Fig. 21 we show the distribu- 
tion of the formation redshifts for model elliptical galaxies. 
We define the formation redshift as the redshift when 50 
per cent (or 80 per cent) of the stars that make up the fi- 
nal elliptical galaxy at redshift zero are already formed. The 
shaded histograms are for model elliptical galaxies with stel- 
lar mass larger than 10 11 Mq, while the open histogram is 
for all galaxies with secure morphology (stellar mass larger 
than 4 x 10 9 Mq). The figure clearly demonstrates that the 
stars in more massive ellipticals are on average older than 
those in their less massive counterparts, but the scatter of 
the distribution is rather large and there is a non-negligible 
fraction of model galaxies whose stars are formed relatively 
late. 

It is important, however, to distinguish the early for- 
mation times of the stars that make up the elliptical galaxy 
population (reflected in the 'down-sizing' scenario) from the 
assembly time of the more massive ellipticals. If massive el- 
lipticals form a large fraction of their stars in a number of 
distinct progenitor systems before they coalesce, these two 
times may well be quite different. Fig. demonstrates that 



2 The actual bin size used for Figs. [I] and is 7 X W x Mq < 
M stars < 2 X 1CP+ 1 Mq, with x = 8, 9, 10, and 11. 

3 The actual bin size used is 7 X W x Mq < M 20 o < 2 X 10^+ 1 Mq , 
with x = 11, 12, 13, and 14. 
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Figure 4. Distribution of the formation redshifts of model ellipti- 
cal galaxies. In the upper (lower) panel, the formation redshift is 
denned as the redshift when 50 per cent (80 per cent) of the stars 
that make up the elliptical galaxy at redshift z = are already 
formed. The shaded histogram is for elliptical galaxies with stel- 
lar mass larger than 10 11 Mq, while the open histogram is for all 
the galaxies with mass larger than 4 X 10 9 Mq. Arrows indicate 
the medians of the distributions, with the thick arrows referring 
to the shaded histograms. Note that more massive ellipticals typ- 
ically form their stars earlier. 



this is indeed the case in our model. We here show the distri- 
bution of the assembly redshifts for the same galaxies that 
we analysed in Fig. 0] We define the assembly time as the 
redshift when 50 per cent (or 80 per cent) of the final stel- 
lar mass is already contained in a single object. For galaxies 
more massive than 10 11 Mq, the median redshift when half 
of the stars are formed is ~ 2.5 (upper panel of Fig.0J, but 
for the same galaxies, half of their stars are typically assem- 
bled in a single object only at redshift ~ 0.8 (upper panel 
of Fig.|SJ. In addition, more massive galaxies assemble later 
than less massive ones, and only about half of the model 
elliptical galaxies have a progenitor with mass at least equal 
to half of their final mass at redshifts > 1.5. The assembly 
history of ellipticals hence parallels the hierarchical growth 
of dark matter halos, in contrast to the formation history of 
the stars themselves. Note that the 'gap' between assembly 
redshifts and formation redshifts for the stars grows towards 
more massive ellipticals. Figs. 2] and imply that a signifi- 
cant fraction of present elliptical galaxies has assembled rel- 
atively recently through purely stellar m ergers. This finding 
agrees with recent observational re sults llvan D okkumll2005t 
iFaber et al1l2005l : iTran et alJl2005l : iBell et al.ll2005l) . 

Table lists, for different mass bins, medians and up- 
per and lower quartiles of the distributions of lookback times 



0.35 l | | | | | 

0.30 r 




redshift 

Figure 5. As in Fig. |I] but for the assembly redshifts of model 
elliptical galaxies. We define the assembly redshift as the time 
when 50 per cent (80 per cent) of the stars that make up the 
galaxy at redshift zero are already assembled in one single object. 
Note that more massive ellipticals typically assemble their stars 
later (cf Fig.gJ. 



corresponding to the formation and assembly redshifts de- 
fined above. 

The large volume of our simulation allows us to study 
how properties of model elliptical galaxies depend on their 
stellar mass and on the environment. Fig. E] shows how the 
luminosity-weighted age (panel a), the metallicity of the 
stellar component (panel b), and the B— V colour (panel 
c), depend on the galaxy stellar mass. In each panel, filled 
circles represent the median of the distributions, while the 
error bars mark the upper and lower quartiles. In the up- 
per panel, the empty circles show the lookback times corre- 
sponding to the formation redshifts as defined in the upper 
panel of Fig. 01 The ages of model elliptical galaxies range 
from ~ 4 Gyr for galaxies with mass a few times larger than 
10 s M to ~ 10 Gyr for galaxies with mass ~ 10 12 M Q . It is 
interesting to note that the lookback time corresponding to 
the redshift when half of the stars had formed is a very 
good approximation to the luminosity-weighted age over 
the full range of masses shown. The age of model elliptical 
galaxies also seems to flatten at stellar masses ~ 10 10 M . 
The same flattening is observed for the B— V colour and, 
although weaker, for the stellar metallicity. Note that the 
scatter in these quantities (particularly for the colour and 
metallicity) is very small, indicating that the main driver of 
these trends is the stellar mass, as also reflected in Fig. 
We note in passing that our model elliptical galaxies follow 
a colour- magnitude relation that is well defined up to z ~ 1. 
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Table 1. Formation times and assembly times for model elliptical galaxies in different mass bins. The first two columns indicate the 
extrema of the mass bins. The next columns list the lower quartile, the median, and upper quartile for each of the formation and assembly 
times defined in the text. Tf50, Tf80 represent the lookback times corresponding to the redshifts when 50 or 80 per cent of the stars 
were first formed. Ta50 and Ta80 represent the lookback times corresponding to the redshifts when 50 or 80 per cent of the mass was 
first assembled in a single object. All times are in Gyr. Masses are in units of Mq. 
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Figure 6. Median luminosity-weighted age (panel a), stellar 
metallicity (panel b), and B— V colour (panel b) of model ellipti- 
cal galaxies as a function of their stellar mass. Symbols indicate 
the median value of the distributions at each mass, while the er- 
ror bars link the upper and lower quartiles. The open symbols in 
the top panel correspond to the lookback time of the upper panel 
in Fig. E| 



The results shown in Fig. |S| also indicate that this relation 
is mainly driven by metallicity, in line with the common in- 
terpretation of the evolution of the slope and the zero-point 
of the observed relation as a function of redshift. We plan to 
investigate the relative contribution of age and metallicity in 
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Figure 7. The same quantities as in Fig. |§] but now calculated 
as mass-weighted average over all galaxies in a halo and shown as 
a function of the virial mass of the halo. Only elliptical galaxies 
with stellar mass larger than 4 X 10 9 Mq are included here. 



shaping the observed colour-magnitude relation in a future 
paper. 

In Fig. Q we show the same quantities as in Fig. but 
as a function of the virial mass of the halo in which the 
galaxies reside. Only galaxies with stellar mass larger than 
4 x 10 9 M are included here. For each halo, we compute the 
mass-weighted average age, metallicity and colour. Fig. |7| 
shows the median of these values and the upper and lower 
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quartiles of the distribution. Elliptical galaxies in high den- 
sity environments are on average older, more metal rich, and 
redder than isolated elliptical galaxies. Elliptical galaxies in 
clusters form a homogeneously old population and ellipti- 
cal galaxies in groups (haloes with mass ~ 10 13 M ) are as 
old as or slightly younger than galaxies in massive clusters 
(haloes with mass ~ 1O 15 M0), while elliptical galaxies in 
smaller groups exhibit a lower median a ge. This is in agree- 
ment with recent observa tional results l|Terlevich fc Forbesl 
120021 ; IProctor et alJl2004T) . The median stellar metallicities 
of galaxies in our clusters are higher than the corresponding 
values for the field. Elliptical galaxies in these systems show 
a remarkably small spread in colour with a median B— V 
colour that is almost independent of environment down to 
the mass-scale ~ 10 13 Mq of groups. 

We can also investigate how the properties of model 
elliptical galaxies depend on cluster-centric distance. This 
is shown in Fig. |H1 based on all elliptical galaxies in haloes 
with M200 ^ 8 x 10 14 Mq and with stellar mass larger than 
4 x 10 9 Mq. We find in total 51 clusters in the whole sim- 
ulation box with virial mass larger than this value. Fig. |H1 
shows that galaxies closer to the centre are on average older 
and more metal rich than galaxies at the outskirts of these 
clusters. The bottom panel of Fig. |H1 shows that this trend 
is partly driven by mass segregation. A radial dependence 
of galaxy properties is, however, also a natural consequence 
of the fact that mixing of the galaxy population is incom- 
plete during cluster assembly. This implies that the cluster- 
centric distance of the galaxies is correlated with the time 
they were accreted o nto a larger structure ijDiaferio et alJ 
l200ll:lGao et al .120041) . The values of the median age, metal- 
licity, and colour all flatten at distance ~ R.200- The scatter 
shown in Fig. 0is rather large but in the very centre of the 
clusters, where elliptical galaxies are all very old (~ 12 Gyr), 
they have about solar metallicity, and have very red B— V 
colour (~ 0.95). 

In the hierarchical galaxy formation scenario, elliptical 
galaxies form through mergers of smaller units, and larger 
systems are expected to be made up by a larger number of 
progenitor galaxies. A very interesting question is therefore 
how large the number of progenitor systems of galaxies is, 
and how this number varies as a function of final mass. To 
get a quantitative handle on this question, we define for each 
galaxy an effective number of stellar progenitors by comput- 
ing the quantity 



iVeff = 



Mi 



2 J^i m i Mif OTm ' 



(3) 



where rrij denotes the masses of all the stars that make up a 
galaxy of final mass Manal = X), mi - The quantity Mi,f orm 
gives the stellar mass of the galaxy within which the star 
i formed, at the time of formation of the star. In the case 
where all stars form in a single object that grows to the final 
stellar mass without experiencing any merger, Eq. can 
be viewed as a discretised form of the integral 



J~2MdM 

which evaluates to N e ff = 1 independent of the detailed star 
formation history. However, if a galaxy is assembled from 
several pieces we expect a larger value of N e g, because then 
the values of Mi f orm that enter the sum in the denomina- 
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Figure 8. Median luminosity-weighted age (panel a), stellar 
metallicity (panel b), B— V colour (panel c), and stellar mass 
(panel d) for model elliptical galaxies in dark matter haloes with 
mass ^ 8x 10 14 Mq as a function of the distance from the cluster 
centre. Symbols and lines have the same meaning as in Fig. [5] 



tor of Eq. © become lower. For example, if the stars of a 
galaxy were formed in two progenitors of equal final size, 
which then merged into a single object without any further 
star formation, we obtain N b h = 2. Note that in more gen- 
eral cases we will obtain fractional values for N e s due to 
the mass-weighting of the progenitors, which is built into 
the definition of N e ff. For example, if a galaxy is made up of 
three pieces that contain one half, one quarter and one quar- 
ter of the final stellar mass, respectively, one gets JV e ff = 8/3, 
which is less than the absolute number of progenitors. This 
reflects the fact that the majority of the stars formed in a 
single object. The mass- weighting hence delivers an effec- 
tive count of the progenitors by giving weight only to those 
progenitor systems that contribute significantly to the final 
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Figure 9. Effective number N e g of progenitors as a function 
of galaxy stellar mass. Filled circles represent the median of the 
distribution in our default model, while the error bars indicate the 
upper and lower quartiles. Empty circles and the corresponding 
error bars are for a model where bulge formation through disk 
instability is switched off. The vertical dashed line corresponds 
to the limit above which our morphological type determination is 
robust (see Sec. I5t. 



stellar mass of the galaxy. In contrast, a count of the to- 
tal number of progenitors would be dominated by the large 
number of negligibly small satellites that fall into a galaxy 
during its growth in a hierarchical universe. We therefore 
argue that N e a is a more useful proxy for the number of sig- 
nificant mergers required to assemble a galaxy. We caution 
however that a straightforward interpretation of N e g in the 
context of spheroid formation is complicated by the fact that 
bulges can also grow in our model without mergers, through 
disk instabilities. 

In Fig.|5] we show the effective number of progenitors as 
a function of galaxy stellar mass. Filled circles represent the 
median of the distributions in our default model while empty 
circles represent the median of the distributions in a model 
where bulge growth through disk instability is switched off. 
Interestingly, bulge growth through disk instability seems to 
be an efficient process for intermediate mass ellipticals but 
rather ineffective for the most massive ellipticals in our sam- 
ple. As expected, more massive galaxies are made up of more 
pieces. Fig. [5]shows that the number of effective progenitors 
is less than 2 up to stellar masses of ~ 10 Mq, indicating 
that the formation of these systems typically involves only 
a small number of major mergers. Only galaxies more mas- 
sive than ~ 10 11 Mq are built up through a larger number of 
mergers, reaching up to N e g ~ 5 for the most massive galax- 
ies. We recall that these most massive elliptical galaxies are, 
however, also the ones with the oldest stellar populations. 
We note that the monolithic collapse scenario would predict 
Neg — 1 for these large ellipticals, in marked difference to 
our hierarchical prediction. 



6 DISCUSSION AND CONCLUSIONS 

We have combined a large high-resolution cosmological TV- 
body simulation with semi-analytic techniques to investi- 
gate the formation and evolution of elliptical galaxies in a hi- 
erarchical merger model. Understanding the formation and 
the evolution of these systems represents an issue of funda- 
mental interest as 50 per cent or more of the stellar mass in 
the loc al Universe appe ars to be in early-type systems and 
bulges (iBell et alJl2003) . 

In this paper, we have focused on the dependence of 
the star formation histories, ages, and metallicities on en- 
vironment and on galaxy stellar mass. In our model, we 
find that elliptical galaxies in denser environments are on 
average older, more metal rich, and redder than the general 
population of 'field' ellipticals. This can be attributed to the 
fact that high density regions form from the highest density 
peaks in the primordial field of density fluctuations, whose 
evolution is somewhat accelerated with respect to regions of 
'average' density. There is also a clear trend for increasing 
ages and metallicities, and for redder colours, with decreas- 
ing cluster-centric distance. This can again be viewed as 
a natural expectation of hierarchical models where the dis- 
tance of the galaxies from the cluster centre is correlated 
with the time they were accreted onto the larger system. 
When this infall happens, we assume that the galaxy is 
stripped of its hot gas reservoir so it is no longer able to 
accrete fresh material for star formation. The galaxy then 
rapidly consumes its cold gas moving towards the red se- 
quence. 

We have also investigated how the properties of model 
elliptical galaxies change as a function of the stellar mass. 
We have shown - and this is perhaps the most important 
result of our study - that in our model the most massive el- 
liptical galaxies have the oldest and most metal rich stellar 
populati ons, in agreement with observational results (see for 
example iNelan et al . 2005). In addition, they are also char- 
acterised by the shortest formation time-scales, in quali- 
tative ag reement with the recently established down-sizing 
scenario JCowie et al.lll996t) . However, these old ages are in 
marked contrast to the late assembly times we find for these 
galaxies. In fact, our results show that massive ellipticals 
are predicted to be assembled later than their lower mass 
counterparts, and that they have a larger effective number 
of progenitor systems. This is a key difference between the 
hierarchical scenario and the traditional monolithic collapse 
picture. 

Our results disagree with previous semi-analytic models 
that found a trend for mor e massive ell i pticals to be younger 
than less massive ones (iKauffmannl Il996bl : iBaueh et alJ 
119961: iKauffmann fc CharlotJll99St) . In order to understand 
the origin of this discrepancy we have re-run our model with 
different assumptions. Fig. llOl shows the average star forma- 
tion histories of all model elliptical galaxies split into bins 
of different stellar mass and normalised to the total mass of 
stars formed, as in Fig. □ In the upper panel, we repeat the 
results obtained for the 'standard' model used in our analy- 
sis, which employs the suppressi on of cooling flows by central 
AGN activity as introduced bv lCroton et alJ ll200Fj) . In the 
middle panel, we show the same results but for a model in 
which no AGN feedback and no artificial cooling cutoff is 
included. Finally, in the bottom panel, we show results ob- 
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tained for a model in which galaxy cooling is switched off in 
haloes with V v ir > 350 kms -1 - the ad hoc s uppression used 
in many previous models, including those of iDe Lucia et ail 
i2004tJ) . 

Fig. 1101 clearly shows that when no suppression of the 
condensation of gas in massive haloes is considered, the most 
massive ellipticals have the most extended star formation 
histories. Too many massive systems are, however, produced 
at redshift zero, at odds with observations. Late mergers 
and late accretion, which still involve a substantial amount 
of gas in this model, cause the formation of luminous and 
young bulge stars. An artificial cutoff of the gas condensa- 
tion, similar to that employed in previous models, produces 
results that are qualitatively similar (see also lDe Lucid 120041) 
to those obtained wi th the more physica lly motivated AGN 
model introduced bv lCroton et al.l i2005l) . as is shown in the 
bottom panel of Fig. QJ3 The fi gure also indicates however, 
that the model does not produce a monotonic behaviour as 
a function of stellar mass. 

This is better seen in Fig. 1111 where we show the me- 
dian luminosity-weighted ages of model elliptical galaxies 
as a function of their stellar mass for the same three mod- 
els. Filled circles show the result for the model with AGN 
feedback, open circles show the result for the model with- 
out AGN feedback and without any artificial cutoff, and 
filled triangles show the result for the model with the artifi- 
cial cooling cutoff. The lines coincide perfectly up to stellar 
masses ~ 5 x 10 10 Mq. For larger masses, the median age 
stays almost constant for the model with AGN feedback, 
decreases for the model without suppression of the cooling 
flows, and shows a non-monotonic behaviour for the model 
with an artificial cooling cutoff. We note, however, that the 
differences between this scheme and a model with AGN feed- 
back are small. In our analysis, a cooling flow cut-off is hence 
able to approximately produce the same result as AGN feed- 
back, which is still different from many earlier results. 

We note that the earlier semi-analytic results were 
based on Monte-Carlo merger trees constructed with the 
extended Press-Schechter formalism, and not by measuring 
merger trees directly from high-resolution numerical simu- 
lations as we have done here. Although there have been sug- 
gestions that the extended Press-Schechter theory may not 
provide a sufficiently accurate description of the merger trees 
feenson. Kamionkowski fc Hassanill2005l) . it seems unlikely 
that this is responsible for reversing the trends found in this 
analysis. We believe this reversal to be a combination of the 
change in the ph ysical model and in the c osmol ogy. We note 
that the work oi'.K.-iiilSiiiaiin fe Charlotl jl99Sfl was carried 
out in the framework of a cosmological model with critical 
matter density, where massive haloes are formed much later 
than in a low-density Univers e normalised to the same clus- 
ter abundance. In addition, iKauffmann fc Charlod (Il998f) 
employed an artificial cooling cutoff corresponding to a crit- 
ical velocity equal to 500 kms -1 . This leaves room for late 
gaseous mergers and gas accretion that can substantially 
rejuvenate the stellar population of elliptical galaxies. 

Another important difference between our model and 
previous ones is that we explicitly follow dark matter sub- 
structures within each halo, even after their progenitor ha- 
los ha ve been accreted by larger structures. ISpringel et alJ 
(2001) have shown that this allows a more faithful tracking 
of the orbits of galaxies and an improved modelling of the 
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Figure 10. Average star formation histories of model elliptical 
galaxies split into bins of different stellar mass, normalised to the 
total mass of stars formed. In the upper panel, results are shown 
for the 'standard' model used for our analysis, which employs the 
suppression of coolin g flows by ce ntral galaxy AGN activity, as 
introduced bv ICroton et all l200fj| . In the middle panel, results 
are shown for a model without AGN feedback and without any 
cooling cutoff. Finally, in the bottom panel, results are shown for 
a model without AGN feedback but with a cooling cutoff with a 
critical velocity equal to 350 kms -1 . Different line styles have the 
same meaning as in Fig. HI 



actual merging rate in any given halo. Most of the galaxies 
we classify as ellipticals are indeed genuine 'satellite' galax- 
ies with their own self-bound dark matter subhalo. If too 
many of these satellite galaxies are assumed to merge on 
too short a timescale, excessively bright and blue central 
galaxies result. 

The short formation time-scales we find for very mas- 
sive ellipti cals are qualitatively in agreement with those re- 
quired by iThomas et al.l <|2005T) in order to reproduce the 
observed a-element enhancements. The detailed census of 
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Figure 11. Median age of model elliptical galaxies as a function 
of galaxy stellar mass. Filled circles show results for our default 
model. Empty circles are for a model without AGN feedback and 
without any cooling cutoff. Filled triangles are for a model with- 
out AGN feedback but with a cooling cutoff with a critical velocity 
equal to 350kms — 1 . 



ages and metallicities for stars in our model elliptical galax- 
ies depends on details of our feedback model and chemical 
enrichment scheme. We plan to come back to these issues in 
future work. 

Our results demonstrate that an apparent 'down-sizing' 
in the formation of ellipticals is not in contradiction with 
the hierarchical paradigm. Modern semi-analytic models of 
galaxy formation do predict 'anti-hierarchical' star forma- 
tion histories for ellipticals in a ACDM universe even though 
the assembly of these galaxies is indeed hierarchical. 
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